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Abstract 

We consider two well-known integrable systems on the plane using the concept of natural Poisson 
bivectors on Riemaninan manifolds. Geometric approach to construction of variables of separation 
and separated relations for the generalized Henon-Heiles system and the generalized system with 
quartic potential is discussed in detail. 



1 Introduction. 

A fundamental requirement for new developments in mechanics is to unravel the geometry that underlies 
different dynamical systems, especially mechanical systems. In fact, geometric analysis of such systems 
reveals what they have in common and indicates the most suitable strategy to obtain and to analyze 
their solutions. 

The Hamilton- Jacobi theory for finite-dimensional Hamiltonian systems is well understood in both 
classical and geometric points of view, see foundational works of Jacobi, Stackel, Levi-Civita and others. 
Apart from its fundamental aspects such as its relation to the action integral and generating functions of 
symplectic maps, the theory is known to be very useful in integrating the Hamilton equations using the 
technique of separation of variables. 

An integrable system is separable if there are variables of separation and n separation relations 



^i{ui,Pui,Hi,. . . ,Hn) =G , i = l,...,n, with det 



dH, 



^0, (1.1) 



connecting single pairs Ui,pui of canonical coordinates with the n Hamiltonians Hi, . . . , H„. Solving 
these relations in terms of Pm one gets the Jacobi equations and a corresponding additively separable 
complete integral of the Hamilton- Jacobi equation 



i—l 



Any separable system is a bi-integrable system [71 [TU], i.e. n functionally independent integrals of 
motion iJj. are in bi-involution 

{H^,Hk} = {H^,Hk}' = 0, i,k = l,...,n, (1.2) 

with respect to compatible Poisson brackets {., .} and {., .}' associated with the Poisson bivectors P and 
P' , such that 

[P, P] = 0, [P, P'] = 0, [F', P'] = 0. (1.3) 

Here [., .] is the Schouten bracket. 

For the given integrable system fixed by kinematic bivector P and a tuple of integrals of motion 
Hi , . . . , Hn bi-Hamiltonian construction of variables of separation consists in a direct solution of the 
equations ()1.2|) and ()1.3|) with respect to an unknown bivector P' [HI [9]. The main problem is that 



geometrically invariant equations ()1.2|1.3p a'priori have infinite number of solutions [3 [10]. So, in order 
to get a search algorithm of effectively computable solutions we have to narrow the search space by using 
some non-invariant additional assumptions. 

The aim of this note is to prove that the concept of natural Poisson bivectors on the Riemannian 
manifolds allows us to properly restrict the search space and to calculate new variables of separation for 
the well-known generalized Henon-Heiles system and the system with quartic potential [21 [21 0] . 



2 Settings 

In this section we recall some necessary facts about natural bi-integrable systems on Riemannian manifolds 
admitting separation of variables in the Hamilton- Jacobi equation. 

Let Q be a n-dimensional Riemannian manifold. Its cotangent bundle T*Q is naturally endowed 
with canonical invertible Poisson bivector P, which has a standard form in fibered coordinates z = 
{qi,...,qn,Pi,.-.,Pn) on T*Q 



P = 



I 

-I 



(2J) 



The Hamilton function for natural system on Q 



H = T + V = ^ gijPiPj + V{qi, ...,(?„) 



(2.2) 



is a sum of the geodesic Hamiltonian T and potential energy V. 

According to |10[ 111] , for the overwhelming majority of known natural Hamiltonian systems second 
bivector P' usually has a natural form, i.e. P' is a sum of the geodesic Poisson bivector P^ and the 
potential Poisson bivector defined by a torsionless (1,1) tensor field A((7i, . . . , g„) on Q associated with 
potential V 

A,. 



P' = P!r 



( 

-A 



V 



J 2 



E 



fc=i 



ki 



dk 



kj 



Pk 



(2.3) 



dqj dqi 

The geodesic Poisson bivector P^ is defined hy n x n matrix H on T*Q and functions x, y and z: 



p^ 



El \ C^n,7c I ^ 9Hife 

^jfcW^^i — y»fc((?)- 



dpi 



-H 



dpj 



Hi 



V 



E 

fc=i 



(9H 



ki 



an 



kj 



dqj dqi 



Zfe(p) 



(2.4) 



In fact, functions x, y and z are completely determined by the matrix H via compatibility conditions 

= [F^,P^] =0. (2.5) 

We can add various integrable potentials V to the given geodesic Hamiltonian T in order to get integrable 
natural Hamiltonians (|2.2|) . In similar manner we can add different compatible potential matrices A to the 
given geodesic matrix H in order to get natural Poisson bivectors P' ()2.3|) compatible with the canonical 
bivector P. 

The definitions of natural Hamilton functions (|2.2I) and natural Poisson bivectors (|2.3p are non- 
invariant because they depend on a choice of coordinate system. Below we show how canonical transfor- 
mations of variables change this definitions of natural Poisson bivectors. 
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3 Generalized Henon-Heiles system 



Let us consider a generalized Henon-Heiles system [3] defined by the following Hamilton function 

Hi = + f <Z2(3g? + 16^2^) + c, {2ql + |) + | + | (3.1) 



and second integral of motion 

TT 4 , /'7?(3cig2 + C2) 4c4 4c5\ 2 Ciq^ cfqf (c2 - 3ci92)(c2 + Ci52)qf 

(3.2) 

, C2C5 + 4c| + 3ciC5(72 , 8C4C5 4c| 
+ C1C492 H 4 ^ — + — ■ 

91 ql qi 

At C4 = C5 = variables of separation have been obtained in [3], see also discussion in [TJ[S]. In this section 
we recover these results in the framework of the bi-Hamiltonian geometry and obtain new variables of 
separation for the generic case. 

3.1 Case C4 = C5 = 

Let us suppose that the required bi- vector P' has a natural form (|2.3p . Substituting Hi 2 (|3.1l3.2p into 
the equation (|1.2[) and solving resulting equation together with (|1.3p one gets two distinct solutions P[ 
and P^. 

First solution P[ is defined by the following geodesic matrix 

H(i) = i I I (3.3) 

and potential matrices 



\P1P2 \pI 



Ad) - I ^ I , (3.4) 

-# cig| + C2gi 

for which 

Xjfc = yjfc = Sjkqk , zfc(p) = 0. 
The corresponing recursion operator A^i — P[P^^ yields the integrals of motion 

^fc-^trATf, fc = l,2, (3.5) 
which form a bi-Hamiltonian hierarchy, i.e. the Lenard relations hold 

P'd'Hi^PdH2. (3.6) 
These integrals are the following functions of initial integrals of motion p.l|3.2p 

Hi = 2i/i 7^2 = + ^ . 

So, the first recursion operator iVi = P-^P^^ gives rise to the action variables Hi^2 P-Sp with trivial 
dynamics. 

Second solution is defined by the matrices 




A(2 



/ CiQa £2 CiQi g2(3ci(72 + 2C2) \ 

'2 4 8 qi ' 

. cigi cig2 . 

^16 4 ^ 



(3.7) 
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and functions 

Xji = Yii = -qi , Zfc(p) = . 
Instead of the Lcnard relations (|3.6p here we have the foUowing relations 

2 



P'dH, = P^F.j dHj, 



where -F is a so-called control matrix 



2 



Ciq2 + C2 

-2 



z = 1,...,2, 



9? 



(3.8) 



16i/2 9i ^ -iplqi ^ + Ciq2 + C2 



Eigenvalues of F coincide with the eigenvalues of recursion operator, which are the desired variables of 
separation. 

So, the second recursion operator N2 = Pl^P^^ yields the variables of separation with non-trivial se- 
parated relations. These variables of separation and the corresponding separated relations are considered 
in the next section. 

As usual, the recursion operators iVi 2 generate two infinite families of solutions of the equations 
(OTT31) 

m = ...,-1,0,1,..., (3.9) 



p(™) _ /\rm p 
^1,2 — ^^1,2 ^ 1 



associated with the Hamiltonians Hi, 2 (|3.1|3.2|) . 
3.2 Case 04,5 ^ 

According to [10] at C5 7^ we have to apply the following canonical transformation 



Pi -^Pi 



-2C5 

1! 



(3.10) 



to the natural Poisson bivectors P[ 3. Namely, in generic case the integrals of motion Hi 2 (13.113.21) are 
in involution with respect to the Poisson brackets associated with the shifted bivectors 



Pi = P[ 



V^5 
if 



I 



pi 

* 



H 

P2 





-2C5 







^ql + Qciql + 4:C2q2 



(3.11) 



and 



27^2^ 



4 



pi + 
* 



1 



-2C5 

if 



P2 





3cig? 







6cigf -I- 4c2g2 




(3.12) 



These bivectors were obtained by substituting non-homogeneous polynomial ansatz for the geodesic bivec- 
tor into the equation (|2.5p . In contrast with the pair of bivectors P{ 2 p.3|3.7p these bivectors Pi. 2 
generate variables of separation with nontrivial separated relations in both cases. 
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Firstly, let us consider the recursion operator N2 = P2P ^ and its eigenvalues ui,2, which are the 
roots of the following polynomial 



B{X) = (A-ui)(A-U2) = A - (^-2 + + -5 -^JA 

(3.13) 

_ ci (4^192 - 2giPiP2 - C2gig2) cf(8g^ + g2) ciV-2c5(4pig2 - giP2) CiC5g2 
16ql 16 8gf 2g8 " 

According to [8l[9l[T0] now we are looking for the auxiliary polynomial A{X) — aiX + ao, which is solution 
of the equation 

{BjX), Aj,)} ^ - ^"^^^ + + - ^'^'^ + + ^(^^ , {A(A),A(,)}.0 (3.14) 

A — /i 

with respect to unknown functions ai^O: '^1,2 and do- In our case ^2 = c?o = 0, do = 1, and desired 
polynomial 

64(g3pi + V-2C5) , 4(4pi(72 - giP2) 16V-2c5g2 

^[■^) — 2T 4 

cfqj aqi ciql 

satisfies the following equations 

{i?(A), A(m)} = (i3(A) - i3(M)) , {A(A), A(/.)} = . 

It entails that 

Pu^ = A{X = uj) , = 5ij , i,j = 1, 2, 

are canonically conjugated momenta. 

An inverse canonical transformation from variables of separation to the initial variables looks like 



4iPu,-Pu2) , „„C2(Ml+U2)-4(u?-UiW2-Mi) 

51 = \ r + 32 5 

y 32 [ui - U2) cf 



Cl(p«i -P«,) \/-2c5 3/ Pui~Pu2 \ , 4(ui+U2)-C2 

64(mi - U2) 91 \32[ui-U2), 



Cl 



V32(ui-U2)/ 4(ui-U2)V 4 
Now it is easy to calculate the separated relations 



$(u,,p„J = <i>+(«,,p.J$_(u,,p„J - '^^'^ ^ ^^^^ + ^'^^P-" = , ft = 1,2, (3.15) 

where 

^±{u,,p^J-^-^-H,±— — + ^- 

These separated relations are given by the afltine equations in Hamiltonians Hi and H2 — iHf . It means 
that the generalized Henon-Heiles system belongs to the Stackel family of separable systems. 

At C4 = C5 = these separated relations p.lSp may be reduced to a pair of distinct separated 
relations 

$+(ui,p„J = 0, and $_(u2,P«2) = • 

and the equations of motion are linearized on two different elliptic curves, see [l]|4j[5]. 

In generic case C4 7^ and C5 7^ equations of motion are linearized on the two copies of non- 
hyperelliptic curve of genus three defined by p.l5p . An explicit description of the linerization procedure 
is an open problem similar to the generalized Kowalevski and Chaplygin systems [51 [3] . 
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Using the first bivector Pi p. lip we can get other variables of separation 

B{X) = ( det(PiP"' - AI) ) ^{X-vi){X- V2) , 
which are related with the previous separated variables by canonical transformation 

cjpl^ 64ul I6C2M?. 

^k = ^^ 2^ + r^, A = 1,2. 

64 cj cj 

4 Generalized system with quartic potential 

Let us consider a generalization of the system with quartic potential [3] with the Hamilton function 

Cl , 4 2 2 , n 4N C2 / 2 , . 2C3 C4 C5 

Hi = — 5 — + ^'^il'i + 892) + y (91 + 492) + -r + 72 + 76 (4-1) 

Z Z H2 HI HI 

and second integral of motion 

4 2/4/i22r.2 4c4 4C5 \ , 

4c^ c^q^ 
+ciq\pl + + 2ciC4<7i + ^ciCAql + Cj^^ + ciC2gi + 2ciC2qlql + 

,2 6 2, 2 4 4, 4ciC3qf /^8c4 4c5 4c2 2ci 12ci \ 

+ci(7i(Z2 + ^9192 + 2 +c5-r + -T2+-;r + -;r + 7472 • (4-2) 

92 V 9i 91 9i 9i 9i92 / 

At C4 = C5 = variables of separation were obtained in see also [5]. In the following section we 
reproduce this result in the framework of bi-Hamiltonian geometry and obtain variables of separation in 
the generic case at 04,5 ^ 0. 

4.1 Case C4 = C5 = 

According to [3 [10], the geodesic matrices 11^^'^^ p.3l3.7p are compatible with another pair of potential 
matrices 



/ ciql , 2 , N 9? , 4 , 2 , C3 ciqlq2 3 , C3 \ \ 

' ^ (3cig2 +C2) + cig2 +C2g2 + — — 7. h 20192 + C2 92 -— 9i 



4 ^ 2 g| 2 V 92' 

ci9i92 4 , 2 , '^a 

-7— Cig2+C292 + -2 

^ 92 



(4.3) 



/'9? , „ 2\ , o2ciqf +C292 - C3 

C2 + Cl — - + 2(72 Ciqi(j2 + 2- 



and 

) V 2 / 9192 . (4.4) 

' Ciqiq2 2 

-C192 

As above, a natural Poisson bivector P[ defined by the matrices 11*^^) and A*^-'^ ^ yields recursion operator 
A^i and integrals of motion "Hi. 2 p.Sp . which form a bi-Hamiltonian hierarchy p. 61) . 

Matrices H^^^ and A'^'^ ' define the second natural bivector the recursion operator N2 = P^^P^^ , 

whose eigenvalues are the variables of separation ui 2 

B{X) = (det(iV2 - AI))'^' = (A-Mi)(A-M2) (4.5) 

^2 ^ 9f ci + 2c2g^ + 2qlciql + 2pl ^ ^ cijqlqlpl + 4c3g^ + 4pfg| - 4pigfgip2) 
29? 492^2 
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In this case equation p.l4|) has another solution c?2 = do = 0, di = 1, and 



AiX) 



cigi 



91 



(4.6) 



It ahows us to calculate the desired momenta 



(4.7) 



In the variables of separation {u,pu) the second integral of motion H2 (|4.2p is a complete square and it 
is easy to prove that Hi and ^/H2 satisfy the following separated relations 



2 ci ci u 



$+(w,p„) ^Hi + -^+2ciupl + ^ + ^^=0, 

2 c\ c\ u 



U = Ml, Pu ^Pui 



U = U2, Pu P«2 J 



(4.8) 



and equations of motion are linearized on a pair of elliptic curves (see PQ SI [S]). 
4.2 Case 04,5 7^ 

At C5 7^ in order to get bi-Hamiltonian structures we apply the same canonical transformation as for 
the Henon-Heiles system p.lOp 



Pi Pi 



-2C5 

9? 



which shifts the Poisson bivectors defined by potential matrices (14.314.41) by the rule 



(4.9) 



Pi=Pi 



1 



-2C5 



« ^^+2\/^ 

* P2 

* * Scigf + (3cig? + 4c2)(72 

* * * 



4C3 

9l 



(4.10) 



and 



P2 =B 



, 2V-2c5 



9? 



+ 



V 



-2C5 

9? 



P2 





8cig| + {Sciql + 4c2)g2 



4C3 

9i 



(4.11) 



At C5 ^ polynomials B{X) and A{X) are obtained from (|4.5I) and (|4.6p using the canonical transformation 
(|4.9p . For instance, a "shifted" polynomial A(A) reads as 



A{X) = 



Pi 



-2C5 

1! 



A qUpi 



1! 



P2q2 



ClQl 
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An inverse canonical transformation from variables of separation to the original variables looks like 



^ 2(mp^^ - U2pIJ ^ 2C3 \ , 2(ui +U2- C2) 



Ul - U2 U1U2 I Cl 



C\{uxVu-, - U2VU2) \/^2c^ / CiUiU2iPui - PU2) C3 



Pi = qi 3 ' 92 



U1-U2 qi y (Mi-M2j '"1M2 

2 f cjcsiuiPu, - U2PU2) , '"1'"2(P«1 / 2^ ^ / n2\\ 

P2 = 7 ^ H ^ 71^ (Ci(p„i -p„J(Mlp„i - U2P«2) - (Ml - M2) 

92 V uiW2("i-'"2) (ui-u2y V yy 

Now it is easy to prove that the corresponding separated relations are equal to 



$(ufe,PuJ = $+(ufc,_p„J$-(ufe,_p„J + {2uk ~ C2)c4 - 2ciV~2c5 UkPu^ =0, fc = 1,2, (4.12) 

where $±(u,p„) are given by (14.81) . The separated relations are afHne equations with respect to the 
Hamilton functions Hi and integral of motion H2 — ^H^. It means that the generalized system with 
quartic potential belongs to the Stackel family of separable systems. 

As above, in the generic case C4 ^ and C5 ^ equations of motion are linearized on the two copies 
of the non-hyperelliptic curve of genus three defined by (|4.12p and we do not know how to solve the 
corresponding Abel-Jacobi equations as yet. 

Using the first bivector Pi (|4.10[) we can get other variables of separation 

1/2 



B{\) = (det(PiP-i- AI)) - {\~vi){X-v2) 



which are related with previous separated variables by the following canonical transformation 

2 ul C2Uk C1C3 

Vk = ciUkP^ ^ H \ , k^l,2. 

Cl Cl Uk 

References 

[1] V. Z. Enolskii, M. Salerno, Lax Representation for two particle dynamics splitting on two tori, J. 
Phys. A: Math. Gen. , v.29, p.L425-L431, 1996. 

[2] B. Grammaticos, B. Dorizzi, A. Ramani, Hamiltonians with high-order integrals and the "weak- 
Painleve" concept , J. Math. Phys., v.25, p.3470, 1984. 

[3] J. Hietarinta, Direct methods for the search of the second invariant, Physics Reports, v. 147, p. 87-154, 
1987. 

[4] V. Ravoson, A. Ramani, B. Grammaticos, Generalized separability for a Hamiltonian with nonsep- 
arable quartic potential , Phys. Lett. A, v. 191, p. 91-95 , 1994. 

[5] S. Rauch-Wojciechowski, A. V. Tsiganov, Quasi-point separation of variables for Henon-Heiles 
system and system with quartic potential, J. Phys. A, v.29, p. 7769-7778, 1996. 

[6] A. V. Tsiganov, Towards a classification of natural integrable systems. Regular and Chaotic Dynam- 
ics, v.ll(3), p.343-362, 2006. 

[7] A. V. Tsiganov, On the two different bi- Hamiltonian structures for the Toda lattice, J. Phys. A: 
Math. Theor. v.40, pp. 6395-6406, 2007. 

[8] A. V. Tsiganov, On the generalized Chaplygin system. Journal of Mathematical Sciences, v. 168, n.8, 
p.901-911, 2010. 

[9] A. V. Tsiganov, New variables of separation for particular case of the Kowalevski top, Regular and 
Chaotic Dynamics, v. 15, n.6, p. 657-667, 2010. 



8 



[10] A. V. Tsiganov, On bi-integrable natural Hamiltonian systems on the Riemannian manifolds^ 
.arXiv:1006.39 14, accepted to Journal of Nonlinear Mathematical Physics, 2010. 

[11] A. V. Tsiganov, On natural Poisson bivectors on the svhere. iarXiv:1010.349 2, 2010. 



9 



